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Abstract 

The general expression for the local matrix t(9) of a quantum chain with the site 
space in any representation of su(3) is obtained. This is made by generalizing t(9) from 
the fundamental representation and imposing the fulfillment of the Yang-Baxter equation. 
Then, a non-homogeneous spin chain combining different representations of su(3) is solved 
by developing a method inspired in the nested Bethe ansatz. The solution for the eigen- 
values of the trace of the monodromy matrix is given as two coupled Bethe equations. A 
conjecture about the solution of a chain with the site states in different representations of 
su(n) is presented. The thermodynamic limit of the ground state is calculated. 



PACS numbers: 75.10Jm, 05.50.+q and 02.20.Sv 



1. Introduction 



The search for integrable spin chains has deserved considerable attention in the last 
years due to the fact that they are interesting physical systems and have a rich mathe- 
matical structure. The best known is the XX Z Heisenberg sit (2) chain with spin S = 1/2 
in every site that gave rise to subsequent development of the quantum groups @-[|J. 
Integrable spin chains with S=l and higher spin chains have been found and solved j5j- 
||ll||.They correspond to higher dimension representations of the quantum group that give 



integrable systems of increasing complexity [12 1-[ 14 



In addition, magnetic hamiltonians can be derived from solution of the Yang-Baxter 
equations (YBE) ( ||15|1|[16|1 ) associated with Lie algebras other than su(2) |I7|] . The solutions 



are found using the Bethe ansatz (BA) for sites with two components or nested Bethe 



ansatz (NBA) for sites with more components [IS]. The introduction of the quantum 



inverse scattering methods (QISM) []19| gave a systematic method to solve those systems. 



The quantum groups give general methods to find new integrable models. 

An interesting problem is to solve integrable chains formed by two kind of states of the 
site. Inhomogeneus solvable models were considered in |[20|| . (see also ||12|| ) . The simplest 
case, an alternating chain with S = 1/2 and S — 1 derived from the su(2) Lie algebra 
was presented in Ref. |H| and in several subsequent works in which the thermodynamic 



properties of these systems was studied p2|-[p5|. 

The system presents interesting features; one of them is that it gives a hamiltonian 
that contains the usual piece coupling pairs of neighboring spins S = 1/2 and S = 1 and 
another piece coupling three neighboring spins. The solution is found using the Bethe 
ansatz. 

In this paper, we are going to solve an alternating chain with the spin of the sites in 
the {3} and {3*} representations of sit(3). We have made an extension of the method used 
in Ref. EJ] for systems where the P and T symmetries are not conserved in order to get 
hamiltonians associated to alternating chains based on the su{2) algebra. 

In a more rigorous sense, we are using the U q (su(3)) algebra and its representations, 
but can be shown that generally for simple algebras g the representations of g and U q (g) 
are isomorphic ||26|| . 

We can obtain two different systems by using as auxiliary spaces the representations 
{3} and {3*}; they will give different hamiltonians, but under a relation between the 
parameters of the local inhomogeneities that we will specified, we can prove that they 
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commute and both systems have the same eigenstates. Then, the more general system will 
be a superposition of those two systems. 

The diagonalization of these hamiltonians requires important modifications of the 



standard method with the NBA [17] [20|. We start building the monodromy matrix in 
the auxiliary space whose elements are operators in the space of states of the chain. The 
main difference is that now we have not a reference state, eigenstate of the operators in 
the diagonal of the monodromy matrix and which is annihilated by all operators under 
the diagonal of this matrix. Then, we introduce a reference subspace in the space of states 
where we can do the second step of the NBA. So, we obtain the equations of the ansatz 
whose properties can be analyzed as in the standard case. The model, since the auxiliary 
space has three dimensions, requires only two steps for the NBA, but the method is easily 
generalizable to more dimensions J27[ . 

The present paper is organized as follows. In the next section we develop the technique 
to obtain the hamiltonians associated to alternating chains ||21|| . In the third section we 
apply the method to alternating chains with the sites in the {3} and {3*} representations 
of su(3). In the fourth section we find the eigenvalues of transfer matrix of the system and 
the equations of our ansatz as a generalization of the NBA. In section five, we analyze of 
the equations of the ansatz and obtain their thermodynamic limit. 



2. Non-homogeneous chain with the site states alternating in two different 
representation spaces 

As is well known, regular solutions of the Yang-Baxter equations (YBE) systematically 



yield integrable chains. In Ref. jzlj an integrable quantum chain with two types of spins 
is described. Following that reference and in order to establish our notation, we are going 
to review how an integrable system follows from a .R-matrix R^(9), which is solution of 
the YBE 

[1 ® R(8 - 0')][R(0) ® l)[l®R(d')] = [R(6')®l][l®R(e)][R(6-6')®l]. (2.1) 
We associate to each site of the chain the t operator 



[taMU d = R b c d a(0), 
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(2.2) 




n a ,b (8)] „, P = 

(c) 



■A 



[ t „,p(9)] a ,b= ^ 

(d) 



P 

b 



Figure 1 



where the indices a and 6 act on the site space and the c and d in an auxiliary space. They 
are shown graphically in fig. 1 (a). Then the YBE can be written in the usual form, 



R{9 - 9') ■ [t{9) <g> t{6')\ = [t(e') <g> t{9)] ■ R(9 - 6'), 



(2.3) 



that graphically is expressed in fig. 2 (a). The <8> product is in the site space and the 
product is in the auxiliary space. 



Equation ( |2.1|) is not the most general YBE. In general we have operators acting on 
pairs of unequal vector spaces. This is represented graphically with lines of different kind. 
We are going to consider two vector spaces denoted by s and a; then we have, besides t, 
the operators t* = R* , i and t represented in fig. 1. They fulfill the YBEs, 



R*{9 - 9') ■ [t*(9) ® t*{B')] = [t*(9') ® t*(9)] -R*(0- 9') 
R*(9 - 9') ■ [t{9) ® t{9')] = [t(9') ® t(9)} -R*(9- 9'), 
R{9 - 9') ■ [t{9) ® t{9')] = [t{9') ® t{9)] ■ R(9 - 9'), 



(2.4a) 
(2.46) 
(2.4c) 
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(a) 



(b) 





(c) 



(d) 



Figure 2 

represented in fig. 2 (b), (c) and (d) respectively. 

In the most general case, we do not require R(9) and R*(6) have P and T symmetry 
nor to be invariant under crossing. Instead, we will assume the following properties: 
i) PT-symmetry, 



ii) Unitarity, 



Ktmc f A-0)=PWa,e6 b ,f, 



iii) Regularity, 



R(0) = c I, 

iv) A matrix M exits such that 

Kp)M b ,eR 9 f %-e - 2rj)Mll oc 6 a , g 8 c , h , 

v) The t-matrices verify, 



(2.5a) 
(2.56) 



(2.6o) 
(2.66) 

(2.7) 
(2.8) 



[ta,b{Q)\otAh,7(-Q)\b,c = p{0)5 ayC 5 an . 
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(2.9) 



We consider a non-homogeneous chain with 2N sites in which the site spaces are 
alternating in the representations {3} and {3*}. This chain has associated the operator 

t^\o, a) = t^eyt^e + a) . . .tg^U^ (*)©. 1>b (* + «) (2.10) 

which is a matrix in the auxiliary space called monodromy matrix, since it describes the 
transportation along the chain. The elements of this matrix are operators on the space 
tensor product of the site spaces. It is graphically represented by fig. 3. 




Figure 3 



Since the t and t matrices fulfill ( |2.3j ), the T^ alt ^ also verifies the YBE 

R(e-e')[T ( - alt \e)(g)T^ lt \e')} = [T( alt \e') ®T( alt \9)}R(9-9'). (2.11) 

Following the standard procedure, we take the transfer matrices 

T^\6,a)=T^\6 1 a). (2.12) 

which are the trace of the monodromy matrices. Due to ( |2.11| ), the operators corresponding 
to different values of the argument 9 do commute, 

[r( alt) (6>,a),T (alt) (6>»] = 0. (2.13) 

The successive derivatives of the transfer matrix at 9 = give us a family of commuting 
operators that describe a solvable system, the hamiltonian of that system being the first 
derivative, 

H=-^\nr^(9,a)\ e=Q . (2.14) 

In a homogeneous chain the hamiltonian is a sum of nearest neighbor interactions 
terms (two-site operators). In our case, it is very different due to inhomogeneities and 
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there are also next-to- nearest neighbor interaction terms (three-site operators). Collecting 
separately the two kinds of terms, the hamiltonian becomes 

2N-1 , 2JV-1 



H _ _L V i _i 1,(2) 

" p» ^ W + c p(a) h h «+W 



add 



3dd 



with 



and 



(\V+l)a,/3;6, 7 = ^a,c(tt)]/3,5[^, 7 (-a)]c,6, 



(\ i+ l )i+2 )a,/3,c;6, 7 ,d = [*a,e («)]/3,5 [*e,d(0)J c,/ fo 7 (-«)] /,&, 

that graphically are expressed in fig. (4. a) and (4.b) respectively. 



(2.15) 

(2.16) 
(2.17) 



i+1 



a a 



( h<1 \ ) = [ t (a) ] [7 (-a)] 

a,P; b,y 1 a,c w J (3,8 1 5,y v 7 J c , b 



(a) 




(2) 



(h i,i+1,i+ 2 U; b , Y ,d ^ V°>W^ W WV^V 



(b) 

Figure 4 



i+1 i+2 



P c 




c 


a C 








* e 




5 


5^aj( 





A similar process can be made by using as auxiliary space the a one. Thus, we define 
the new monodromy matrix 

f<$\e, ^ = #) ai (6> + *K^ 2 (e) . . 4^:^(0 + °V£?\A e )> (2-i8) 

graphically represented in fig. 5. It fulfills the YBE 

R* (6 - 6') [f^ (6-<t)® T (alt) (9' - a)] = [f^ (6' -a)® (6 - a)]R*(6 - 9') . (2.19) 
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T (al P a \ 
T a/ 9 ' ) = 




Figure 5 



The hamiltonian obtained from this monodromy matrix using a formula similar to 
(Q)is 



with 



H 



. 2N 1 27V 

1_ V I 1 V h {2) 



i=2 
i — even 



i = 2 
— even 



(2.20) 



(2.21) 



and 



- /ry\ .. . . _ 

(/l iji+ l >i+ 2)o,o,/i;/9,6,i/ = [*a,« C 47 )]a,cK, i ,(0)] M ,p[t c ,6(-cr)] Pi/3 . 



(2.22) 



The monodromy matrices T^ 1 ^ and T^ 1 ^ fulfill the following YBE 

[i a , b (o -0' + i)] a Xf\^i)T^\e\ -7) = fW(9', -i)T^\d, 7 )[t djC (0 - 0' + 7 )] M , 

(2.23) 

that graphically is expressed in fig. 6. 



e-e'+y 




Figure 6 
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As a consequence of (|2.23 ) the transfer matrices r^ alt ^ and r^ alt ^ commute 

[r (alt) (9, a) , r ( alt ) (0', -a)} = 0, (2.24) 



and then, the derived hamiltonians H and H also commute. Thus, both can be simulta- 
neously diagonalized with common eigenstates. 



3. Quantum chain with the site states alternating in two representations of 

su(3) 



We describe in this section a non-homogeneous chain that we form alternating two 
representations of su(3). We denote a representation by the indices of its associated Dynkin 
diagram (7711,7712). The vector space s is taken as the representation (1,0) = {3} and the 
space a is the generic representation (7ni,7n 2 ). 

The t operator acting on the site and auxiliary spaces, both s, [^7[]can be written [^7 



fi(X 3 q- N 



t(*,7) 



V A 



ei,e 2 



\- 3 q NP ) 



?2 



Or 1 -*?) 



A- 



L /2 



where the parameters A and q have been taken as the functions of 9 and 7 



and the iV matrices are 



(3.1) 



(3.2) 



N a 



2 , k ! r 
3*1 + 3^+3*. 

! L U ! r 
— /ii + -hi + -I, 

3 3 3 



iV 7 = hy - -h 2 + -I, 

3 3 3 



(3.3a) 
(3.36) 
(3.3c) 



where {e^, /i, Q ±/li }, i = 1,2, are the Cartan generators of the deformed algebra U q (sl(3)). 
To obtain the operators i(9, 7), we take ( |3.1[ ) as a basis and write 



i(A) 



A _l (g^-g) 



A (£li Z £) Fi 



A" 



V 



±(\ 3 q- N -\- 3 q NP ) 



l (g F, 
2 ^3 



\ 



AT" 1 



(3.4) 



A-V 7 )/ 



where the operators {Ei,Fi}, i = 1,3, are unknown and will be determined by imposing 
the YBE, 

R{6 - 6', 7) • [t(e, 7) ® W, 7)] = 7) ® 7)] • - <?', 7) , (3-5) 



i?(A, //) 





[to 


,b{0,l)\c,d 


is £ 


jiven 


by 


/a 























°\ 





d 





6 























c 











b 











6 





c 





























a 





























d 





b 











6 











d 























b 





c 





Vo 























a) 



with 



a(A, n) 

c(X,fi) 
d(X,fi) 

The relations obtained are 



(AVV A -A-V«), 

(aV 3 -a-V), 
(ff- 1 -^- 1 , 

(« _1 -«)A-V- 





= q~ 1 q Na E 1 , 


£^ 


= qq N0 E u 




= qq Na F u 




= q~ 1 q N0 F 1 , 




= qq N "E 2 , 


E 2 q N * 


= q- X q N *E 2 , 


F 2 q N " 


- n~ 1 n N " TP 

— q q *2, 


F 2 q N0 


- nn N>3 TP 

— qq f2, 


[Ei,F 1 ] 


= (q' 1 -?) (q 


[E 2 , F 2 ] 


= (q' 1 -q) (q 



—N° 



N a -N fJ 



)■ 
)' 



E3 

F 3 



(q- 1 - q) 
(q- 1 - q) 



q- N0 [E u E 2 ], 

<r [f 2 , Fi], 



(3.6) 



(3.7a) 
(3.76) 
(3.7c) 
(3.7d) 

(3.8a) 
(3.86) 
(3.8c) 
(3.8d) 
(3.8e) 
(3-8/) 
(3.8o) 
(3.8a) 
(3.8i) 

(3-8j) 
(3.8fc) 
(3.8/) 
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and besides, the modified Serre relations 



q- 1 E 1 E x E 2 ~{q + q~ 1 )E 1 E 2 E 1 + qE 2 E 1 E 1 = 0, (3.9a) 

qE 2 E 2 E 1 - (q + q- x )E 2 E x E 2 + q~ l E x E 2 E 2 = 0, (3.96) 

q- x F l F l F 2 - {q + q- x )F x F 2 F x + qF 2 F x F x = 0, (3.9c) 

qF 2 F 2 F x - (q + q- 1 )F 2 F 1 F 2 + q- 1 F 1 F 2 F 2 = 0, (3.9c?) 



should be verified. It must be noted that that the relations ( |3.8a Z| ) are the usual ones 
for the quantum group U q (sl(3)) while the relations ( |3.9a-d| ) are not the usual ones for 
the said group, and because of this, the YBE is not verified if the generators and 
pertaining to the deformed algebra, are taken as Ei and Fi. This induces us to take 

F * = l( ( l~ 1 ~Q)Zifi, (3.10a) 
E i = ^{q- 1 -q)e i Zr\ i = 1, 2 (3.106) 

where and fi, i = 1,2, are the generators of U q (sl(3)) in the representation (mi, m 2 ) and 
Zi are two diagonal operators that were obtain by imposing the verification of the relations 
( p.8a-/|) and ( |3.9a-cj ) . In this way, one obtains the general form of these operators given 
by 

Zi = g«i^-|^+a 3 I ) (3.11a) 

Z 2 = q i h i+(^i+^)h2+b 3 i _ (3.116) 

The knowledge of the operator i permits us to build the monodromy operator of any 
multistate chain that mixes two representations. As an example, for the chain that mixes 
the {3} and the (mi,m 2 ) representations the monodromy operator is 

that is represented graphically as shown in fig. 3. 
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4. Bethe ansatz equations of the models with space sites in different represen- 
tations of su(3) 



In this section, we are going to solve an alternating chain that mixes the {3} and 
{3*} representations of su(3) and the results will be generalized to chains that mix two 
arbitrary representations. 



In this case, the t operator is given by (3.1) that can be written in matrix form 



t(8) 



fa 























°\ 





b 





c 























b 











d 











d 





b 





























a 





























b 





c 











c 











b 























d 





b 





Vo 























a) 



(4.1) 



with 



a(#)=sinh(-0 + 7 ), 

&(#) = sinh(^), 

c(6>) = sinh (7)e 5 , 
d(Q) = sinh (7)e~2~. 



(4.2a) 

(4.26) 
(4.2c) 
(4.2d) 



In the same way, the t is obtained from ( |3.4| ) by taking in ( pMUa, b ) the generators of 
su(3) in the {3*} representation, 





















-1 












Besides, we must fix in ( p.lla, Zj ) the values of a±, as and 63. By taking 



(4.3) 



and rescaling 9 by 



a\ = g, a 3 = 0, 63 = 0, 



= 9+ -7, 

9 1 



(4.4) 



(4.5) 
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we find 



with 



t(0) 



/a 











c 











d \ 





6 





























b 





























6 

















d 











a 











c 

















6 





























b 





























b 





Vc 











d 











a) 



(4.6) 



a(0)=sinh(^+!), 

6(0) = smh (|(6H- 7 )), 
c(6>) = — sinh yy)e 2 



= — sinh ('j)e 2 



(g+r) 



(4.7a) 

(4-76) 

(4.7c) 
(4.7d) 



As we take the {3*} representation as auxiliary space, the i? matrix is 



fa 0000000 0\ 
0ci0600000 
00c000600 
060c00000 
R(0) = OOOOaOOOO 
00000d060 
006000^00 
0000060c0 
Vo 0000000 a) 



(4.8) 



and in this case the matrix M defined in (2^8) is the identity and the functions defined in 
the equations ( |2.6a, 6j) and (|2.7|) are 



p{6) = p*(8) = sinh ( 7 - h) sinh (7 + |(9) s 

c o = Cq = sinh 7, 

p{0) = - (cosh (3 7 ) - cosh (36*)) . 



(4.9a) 
(4.96) 
(4.9c) 



We group two neighbor sites in the chain and form the operator 



odd. 



(4.10) 
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The monodromy matrix that correspond to this model 

T { S ] (0,a) = fg# (0, a)i%% (9,a)... t^'^ (*,«)• 
This operator can be written in the auxiliary space as a matrix 



(4.11) 



A(0,a) B 2 (9,a) B 3 {9,a) 
T alt (6,a) = | C 2 (9,a) D 2 , 2 (6,a) D 2 , 3 (9,a) 
C 3 (9,a) D 3 , 2 (9,a) D 3 , 3 (9,a) 

whose elements are operators in the tensorial product of the site spaces, 

i=odd 



(4.12) 



(4.13) 



s iji+ i being the tensorial product of site spaces (i) and (i + 1) and isomorphic to the {3} 
and {3*} representation product . 



Si,i+i = Si <S> s i+ i ~ {3} <S> {3*}. 
The YBE for can be written in terms of its components 



(4.14) 



B{9) <g> B{6') = R {2) {9 - 6') ■ {B{9') <g> B{0)) = (B(9') <g> B{9)) ■ R {2) {9 - 9'), (4.15a) 
A(9)B(9') = g(9' - 9)B(9')A(9) - B(9)A(9') ■ f (2) (6>' - 9), (4.156) 
D{9) ® B(9') = g(9 - 9')(B(v) <g> D{9)) ■ R^(9 - 9') - B(9) <g> (r (2 )(6> - 9') ■ D{9')), 

(4.15c) 



where 



RW{9) 



/I 








o\ 





d 

a 


6 

a 








b 
a 


c 
a 













1/ 



.(2) 



/i_ 
h+ 



r (2 \9) 



h+ 
h- 



and 



a(6>) 
W) 



h+(9) 



c{9) 
b(9) 



h-(0) 



d{9) 



For the site states, we use the notation 



u 







d 



u 







d 








(4.16) 
(4.17) 



| . (4.18) 
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In order to find the eigenvectors and eigenvalues of 



we find inspiration us in the NBA method and look for an eigenstate of A that serves as a 
pseudovacuum. For this purpose, we build the subspace of sui + x) generated by the vectors 
| u, s > and | u,d>, that we call lOj, and then the subspace 

SI = W! <g> ws <g> • • • ® w N (4.20) 

of the total space of states of a chain with 2N sites. 

In a non- homogeneous chain, we have not a state || v > such that 

D itj || v >oc 5ij || v > . (4.21) 

For this reason, the NBA method cannot be used. Our method, instead, starts with a 
state || 1 >G fl verifying 

A{6) || 1 >= [a(d)] N *$(e)] N Z || 1 >, (4.22a) 
Bi||l>^0, i = 2,3, (4.226) 
d\\l>=0, i = 2,3, (4.22c) 
Dij || 1 >G O, i,j = 2,3, (4.22c?) 

N3 (N£) being the number of sites in the representation {3} ({3*}). In order to simplify 
the exposition of our method, we take ./V3 = iVjf = N. 

Following the steps inspired in the NBA, we apply r-times the B operators to || 1 > 
and build the state 



=^(in,---,ii r ) = B^dn) ■ ■ ■ B ir (/i r )X iu ... jir || 1 >=B(jj, 1 )®...®B(jjL r )X || 1 >, 

(4.23) 

Xi lr .. ) i r being a r-tensor that, together with the values of the spectral parameters 
H\, ■ ■ ■ , fj, r , will be determined at the end. 

The action of A(n) and Di^n) on \1/ is found by pushing them to the right through 
the Bi j (/Xj)'s using the commutations rules ( |4.156, cj ). Two types of terms arise when A 
and Di j pass through S's: the wanted and unwanted terms, similar to obtained in the 
NBA method. The first one comes from the first terms of ( [4.156, cj )- In this type of terms 
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the A or Da and the S's keep their original arguments and give a state proportional to 
The terms coming from the second terms in ( [4.155, cj) are called unwanted since they 
contain Bi((i) and so they never give a state proportional to so, they must cancel each 
other out when we sum the trace of T alt . The wanted term obtained by application of A 
is 

r 
3=1 

and the k-th unwanted term 

r 

-[a(n k )] N3 [b{n k )] N z \[g{n3-Hk) (fl(u)f< 2) (/z fc - u)) ® B(ji k+1 ) <g> • • • 

• ■ • ® S(// r ) ® ® S(^ fc _ 1 )M<*- 1 )X || 1 >, 

(4.25) 

M being the operator arising by repeated application of ( |4.15a[ ) , 



B(m) ® • • • ® S(/x r ) = BGufe+i) ® • • ■ fl(A*r) ® B(a*i) • • • ® B^-ijM^" 1 ). (4.26) 

The application of the operators Dij(p) to the state \I/(/7) is a little more laborious 
but straightforward. The wanted term results 

r 

[AsjM-BiiM • • •■B ir (/i r )X il) ... >ip || 1 >] wanted = JJfif(// - n^Bj^x) ■ ■ ■ 

■ ■■B J MR?j;r> ir (v - ^) • • • Rfj? 2 ' i2 (» - ^) • flg^o* - /'i)/^,-v„.....,, ii i >, 

(4.27) 

where the i?( 2 )'s product is taken in the auxiliary space and has the form 

Qfafla ^Jlf^'^.-flf^^Jif^^ S M '5V (4-28) 

The action of -Dfcj with fc 7^ j on || 1 > is not zero. This is the main difference with 
the models that can be solved by NBA. Then, we try to diagonalize the matrix product 

FM = DM-*M=(ffi ( ^) (4 ' 29) 

By taking the terms in ( |4.27| ) with k = j and adding them for k = 2 and 3, we obtain the 
wanted term 

r 

]\9(P- N) B ii(j i l) , '' B irM^2)(l i ifl X ii,"-,ir II 1 >' ( 43 °) 

3 = 1 
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where 

t (2) (m,/2) =trace(F) = jl) + fl). (4.31) 
In the same form, the k-th unwanted term results 

r 

- Y[ d(Vk - fJ>j) [B(fj 1 )r {2 \fj 1 - fj, k )j <g> B((j, k +i) ® • • ■ 

V J (4.32) 



• • • ® S(/i r ) ® ® B^ k _ x )M {k l) T {2) (iJL k , jt)X || 1 > . 

The sum of the wanted terms and the cancellation of the unwanted terms give us the 
relations 

r ( a)(/i, || 1 >= A( 2 )(/x, /x)X || 1 > (4.33) 

and 

A (2) ( Mfc ,/2) = [aW]^[6( Mfc )]^3 f] (4.34) 

We must now diagonalize ( |4.33|) . 

The state || 1 >G O and the tensor Xi l! ...j r , (ij = 2,3) lies in a space with 2 r di- 
mensions, tensorial product of r two-dimensional spaces Ci,l = 1 • • -r, generated by the 
vectors 

' ' ( ^ l = l---r. (4.35) 



u / 1 \ / 1 

Then, the vector X || 1 > yields in a space with 2 r+N dimensions. In this space, we 
take the element 

|| 1 > (2) = e\ ® e\ ■ ■ ■ ® e* ® \us >i ® • • • ® \us > N , (4.36) 

which is annihilated by C^ 2 \n, jT). (Note that the operators a, (3, 7 and 5 of (|4.28[) act on 
the first part of || 1 >^ and the operators Dij on the second part ). The application of 
the operators and gives 

A^^fl) II 1 > (2) = U»)] N3 [b(ri} N; II 1 > (2) , (4.37a) 
I^W) || 1 >^=n , 1 M »)] N *m] N; II 1 > (2) > (4-376) 

The important fact is that -F(/i, /7) verifies the YBE with the it^ 2 ** matrix given in ( ^.16|) , 

RW(H- //)[F(/i,/2)®F(//,/2)] = [F{n',p) ®F(ii,{l)}R i2) (v-v'), (4-38) 
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which, in a second step, permits us to solve the system. From this equation, we obtain the 
commutation rules 



In this second step, we build the vector 



^ 2 \\fl) = B^(X U fl)---B( 2 \X s ,fI) || 1> 



(2) 



(4.40) 



The action of A^(X,jl) on > gives the wanted term 

[b(X)] N ^X)] N ^f[g(X l -X)B^\X 1 ,fl)...B^\X s ^) || 1 >( 2 \ 

i=l 

and the k-th unwanted term 

s 

-h+(X k - \)[b(\ k )] N > [b(X k )] N S n g(Xi - A fc )5 (2) (A, ^)S (2) (A fc+ i, P) ■■■ 

i=l 
i^tk 

. . . S (2) (A fc _i, /I) || 1 > (2) . 
In the same form, the action of D^(X, jl) on \I/( 2 ) gives the wanted term 

[6(A)] iV3 [a(A)] JV 3 f[g(X - A,) f[ * i? (2) (Ai, pi)... B^(X S , fl) || 1 >^ ., 
i=i j=i ^ 

and the k-th unwanted term 



(4.41) 



(4.42) 



(4.43) 



-h-(\ - X k )[b(X k )] N *ia(X k )] N S f[ g(X k - A,) n 7T l ~ . B<- 2 \\,j!)BW(\ k+1 ,j!) . . . 

i=l ? = 1 3\ A k V>j) 

...s( 2 )(A fc _!,/2) ||1>^. 

(4.44) 

The cancellation of the unwanted terms and the sum of the wanted terms give us the 
equations 



q(Afc) 
b(X k ) 



N* r 



n 



i 



n 



g(Xj - Afc) 

y ^(Afc - m) f = i #(A fc - A,) 
17 



fe — 1, . . . , s, 



(4.45) 



and 



i=l 



Then, by comparing equations (|4.34[ ) and ( |4.46[ ) and calling g{6) 
obtain the coupled Bethe equations 



i=i 



8 



(4.46) 
a(6>)/&(6>), we 

(4.47a) 
(4.476) 



3*1 



and the eigenvalue of the trace of T^ alt ^ 



i=i 



f r~ S ST" 

3 II 9(M - N) \{9{K + M»)] N3 ll9(» - A,) n TTTT^ 

7 = 1 H'j). 



i=l 



i=l 



(4.48) 

that is the solution to the spectrum of our problem. 

The hamiltonian of the alternating chain can be obtained with ( |2.16| ) and ( |2.17| ) . The 
results for and are 



h 



(i) 



sinh 7(1 + 2 cosh 7) 
2(cosh(37) - 



(4.49) 



a=l 



and 



a 

i+2 



a=l 



a=l 



+ fc (A? <8> I i+ i <g> A* +2 - Af ® J <+ i ® Af +2 ) + k'f iA+hi+2 , (4.50) 

where we have used the Gell-Mann matrices A and A for the {3} and {3*} representations 
respectively, being the coefficients, 



sinh 7(1+2 cosh 7) - f / o 

2(cosh(3 7 )-l) 11 " 1 °> 

sinh7(-l+4cosh 2 7) - f _ o 

2(cosh(3 7 )-l) 11 " ~~ °' 



(4.51a) 
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{sinh 2 J (1+2 cosh 7) (3+2 cosh 7) .„ / q g 

2sin h7(cosh(37)-l) 11 tt r <i ) ! 'i ,s 
sinh 2 ^(1+2 cosh 7)(3+cosh7+cosh(2 7 )) . f _qo (4.010J 

2sin h 7 (cosh(37)-l) 11 a ~ d ' °' 

\/3sinh 2 ^(1 + 2 cosh 7 ) 2 
* = 4(cosh(3 7 )-l) ■ (451C) 

sinh7(cosh(37j — 1) 



The term / iji+ i ii+2 is 



v i+2 



+ E K« (Af ® A? +1 ® A? +2 - A? ® A? +1 ® Af +2 ) 

+ w 8 , a (Af ® A? +1 ® A" +2 - Af ® A? +1 ® Af +2 ) 
+f 3)a A« ® Af +1 ® Af +2 + u 8>a A? ® Af +1 ® Af +2 } + z (Af ® Af +1 ® Af 
+Af ® Af +1 <g) Af +2 + Af ® Af +1 ® Af +2 - Af ® Af +1 ® Af +2 ) , (4.52) 



v i+2 



where 



sinh 7 . , 

«J 3 =— ^(2 2 -1 -1 -1 -1 0), (4.53a) 

sinh 7 , , 

w 8 =- — -(0 -1 -1 1 1 0), (4.536) 

— sinh^ — 

v 3 = ^(0 -1 -1 1 1 0), (4.53c) 

2 

v 8 = ^(2 2 -1 -1 -1 -1 0), (4.53d) 

2v 3 



sinh 2 7 



(4.53e) 



v/3 ' 

dfi,u, P are the totally symmetric structure constants of ££7(3), and e M)i , )P is the totally 
antisymmetric tensor. 

As a first generalization, we can apply now the the method to a chain that mixes 
the (1,0) = {3} and (mi,m 2 ) representations. In this model we take again the (1,0) 
representation as auxiliary space; then we have the same i?-matrix ( |4.8|) for the YBE. 

The highest weight of the (mi, m 2 ) representation is 

2m 1 + m 2 m 1 + 2m 2 
A?j = ax H a 2 , (4.54) 
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where a.\ and a 2 are the simple roots of su(3). 

Through ( |3.3a-c| ), (|3.4|), together with the commutation rules of su(3), it is possible 
to know the action of elements of the t matrix on the highest weight vector. We obtain 

t 1 ,i( e )\ A h> = a(0)\A h >, (4.55a) 

*2,2W|A h > = 6~i(0)|A h >, (4.556) 

F 3 ,3(0)|A fc > = &2(0)|A fc >, (4.55c) 

where 

3 2 11 

a(0) = sinh(-6> + (-mi + -m 2 + -) f y), (4.56a) 

h(d) = sinh(-0 + (--mi + -m 2 + -) 7 ), (4.566) 

62(6*) = sinh(-0 + (—mi - 3^2 + 3)7)- (4.56c) 

As before, we group neighbor sites and build the monodromy operator T that can be 
represented by a matrix in the auxiliary space as in ( [1.12|) . The two sites space is now 

s M+ i ~ (1,0)® (mi,m 2 ). (4.57) 

In this space, the subspace Wi is now generated by the highest weight of the (1,0) repre- 
sentation and the subspace V generated by the states 

{|A fc >,/ 2 |A fc >,/f|A fc >,...}, (4.58) 

where / 2 is the generator of sl(3) in the (mi, 7712) representation. 

We form the subspace Q as in (|4.20|) and built the state || 1 >G O which must satisfy 

A(9) || 1 >oc|| 1 >, (4.59a) 

D i}i (0) || 1 >oc|| 1 >, z = 2,3, (4.596) 

Bi \\1>^0, 7 = 2,3, (4.59c) 

d || 1 >= 0, i = 2, 3, (4.59d) 

D iij \\l>en, h 3=2,3, i^j. (4.59e) 

Then, the states ^(fl) analogous to (|4.23|) are 



*(/2) = *(//!,•• -,/i r ) = B^fai) ■ ■ ■ B ir (nr)X ilt ... iir II 1 >=B(in)®...®B(ii r )X || 1 > . 

(4.60) 
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As the YBE depends on the R matrix, we have for the new monodromy matrix the 
same commutations rules ( 4.15a-cj ) as before; then we can repeat the same steps, the only 



difference being in the action of the operators of the monodromy matrix on the state || 1 >. 
The BEs that we obtain in this case are 

[gMfM^r = n #3^4 - ^ ( 4 - eia ) 

j=1 9 Km A*fcJ i=1 

[HW\ N = II 9(Xk - N )f[ 9i ( X :~\ k l (4-616) 
f=i 7=7 9 {M - M) 

where fa, i = 1, ■ ■ • , r, and Xj, j = 1, • • • , s, are the roots of the ansatz, the function g is 
given in ( |4.17f) , and 

U*> = gg, (4.62a) 
*W = g. (4.62.) 

The procedure can be generalized to chains that mix non- fundamental representations 
(mi, m^) and (m^m^), irrespective of the number of sites and their distribution in the 
representations. For this purpose, it is necessary to build the monodromy matrix follow- 
ing an analogous process to used before. If we use a dashed line for the representation 
[vn! \,vn! %), the monodromy matrix T gen (#) can be represented graphically as shown in 
figure 7. 



1 2 




b 



Figure 7 



By calling N± and N2 the number of sites in the representations (mi, m<i) and 
(rn'i, m! 2) respectively, we find the BE for this general chain 

foMl^fcO*)]* = f[ g ^ k ~^\ f[9(K-^) (4.63a) 

j^k 

M*k)\ Ni M*k)] Na = n 9(\ k - Mi) n 9 n " \i ( 4 - 636 ) 

jA 7=1 f ( A * - A *) 
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where g\ and g^ are given in ( [4.62a, and gi and g 2 are the same as the previous ones 
with (mi, 7712) replaced by (m[,m 2 ). 

In the light of this, the generalization for the case of mixed chains with more than two 
different representations seems simple, although the physical models that they represent 
will be less local and the interaction more complex. 

Also we can conjecture about the solution of a non-homogeneous chain combining 
different representations of su(n), each representation introduces (n— 1) functions similar 
to ( [4.62a, fc|) (that we call source functions). The BE are obtained applying the MBA with 
(n — 1) steps, then each solution will have a set of (n — 1) equations (the same number of 
dots in its Dynkin diagram). The first member of the equations will be a product of the 
respective source functions powered to the number of sites of each representation and the 
second a product of g functions coming from the YBE similar to ( [4.63a, 6|) . 



5. Thermodynamic limits of solutions and analysis of Be-the equations 

In this section, we are going to discuss the solutions of the {3}-{3*} model given by 
the equations ( |4.46| ) in the limit for very large N. For that discussion, it is convenient to 
set the parametrization of the spectral parameters 

-fij = iv) - -, (5.1a) 

! Ai = ^f -7, (5.16) 

and N = N 3 + N£ the length of the chain. 

Using such parametrization, the Bethe equations ( [4.47a, fc| ) can be written 



• / (2) 

sin [v k 


+q) 


■ l (2) 

sin [v k 


-q). 


■ i (i) 
sin [Vf. 


-qi 


■ t (i) 
sm(vl 


+<?). 



-in; 



■ I (2) 

sin [v\ 



sm 



n 

r 

n 

j=i sin (v y k - v) 



H\ s • ( (2) (2) 
l 2> TT sm \ v i v k 



sin iv 



(2) 
k 

(1) 
k 



V 



( v k' - v j + q) " sin ( v ^' - v k' + *t) 
(i) 



n 



17) 



,(2) 



,( 2 ) 



(5.2a) 



17) 



n 



• / (2) (1) -i- 

sin(v> ' -v k -ij, 



+ 17) »=i sin ( v i 



+q: 



, (5.26) 



In this regime, the roots must be considered in the interval (— 7r/2, tt/2). Then, we 
define the function 

s sin(x + ia) . . 

<P{X,a) =im-^- f (5.3) 

sm (x — iol) 
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and taking logarithms in (|5.2a, we obtain 



,(2) 7- 



3=1 
r 



,(2) .,(1) 7< 



J2Hv { k ] - = 27r/f \ 1 < k < s(5Aa) 



2) 



(5.46) 



i=l 



where and 7 J [ 2 ' ) are half-integers. 

In the thermodynamic limit iV — > oo, the roots tend to have continuous distributions. 
Unlike what happens in other cases, we cannot distinguish between the roots coming 
from the different types of representations, this can be noted by simple inspection of the 
equations of the ansatz. Due to that, we define two root densities, one for each level, 

1 



pi{vf>) = lim 



iV 3 -oo ^(0 _ (0) 



J = 1,2, 



(5.5) 



Let it be 

Z Na (v) 



1 

Zn* (v) = — 



(i) 



1 (2) 7 

- jj. , — 

jV 3 f-7 TV J '2'J 



AT, 



(2) 



1 



3 J= i 



7) + iV^S 



3 i=1 



(5.6a) 
(5.66) 



The no-holes hypothesis for the fundamental state establishes 



r(i) _ r(i) 
1 k-l 1 k 



i = l,2, 



for all k, 



(5.7) 



that implies 



ZnM 1} ) 



j(D 

iV 3 : 

(2) 



(5.8o) 
(5.86) 



In the thermodynamic limit and for the fundamental state, the derivative of these 
functions are 



AT 
7W 



P2(u). 



(5.9a) 
(5.96) 
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Using the approximation 



1 



N 3 ^qo N3 



Vj ') ~ 



d\f(X)p k (X), 



(5.10) 



together with (|5.9a, |) and ( |5.6a, 6j ), we obtain the system of equations 



Pi(A) 



2tt 



i'(A - n^)p x {n)d^+ 



+ 



2 ^ .-y 



(5.11a) 



02(A) 



2tt 



A^ (A '2' 



0'(A - fi,'y)p 2 (n)dfj,+ 



+ 



that can be solved by doing the Fourier transform, 

00 ^ 

(p(\,a) = n + 2A - % ^ — 

n= — 

Pi (A) = £ ^^^^ 



,2imA — 1\ra\a 



m = — 00 



rn 



[m] 



(5.116) 



(5.12a) 



(5.126) 



Introducing these expressions in the integral equations ( |5.11a, 6|) , we obtain the densities 
in the Fourier space 



pi(m) 
p 2 (m) 



N3 sinh (27|m|) AT| sinh(7|m|) 



when m/0, and 



N sinh(37|m|) A" sinh(37|m|)' 
sinh(7|m|) N£ sinh (27|m|) 
A" sinh(37|m|) A" sinh(37|m|)' 

2(2Ar 3 + AT 3 *) 



Pi(0) 
P 2 (0) 



3N 
2(2A^ 3 * + A^ 
3N 



(5.13a) 
(5.136) 

(5.14a) 
(5.146) 



for m = 0. We note that for N£ = we have again the known result for a homogeneous 
chain. It is interesting to notice the complementarity of the solution for N£ = and the 
solution for ^3=0, 



h( m ) \n 3 =o= h( m ) \n*=o, 
P2(m) |tv 3 =o= Pi(m) \n*=o ■ 



(5.15a) 
(5.156) 
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In the case A3 = N£ = N/2, that corresponds to our alternating chain, the densities 
are given by 

cosh ^7717 

pi{m) = p 2 {m) = — — I . (5.16) 



The free energy is defined by the expression 



lim /(0, 7 ) = -ll g A(0). (5.17) 

N— »oo iV 



Then, taking the dominant term in A(#), 



a+{B) = [a(e)] N *[l(e)] NS ( 5 - 18 ) 

3=1 



the energy is given in this limit by 



m 7) = ~ lg (a(0)) - ^ lg W)) + ^ E + ^' 2 } - (5 ' 19) 

i=i 

Doing the change of variable w = 36>/2 and using equations ( |5.10| ) and ( |5.13a, jp , the 
free energy can be written in the more transparent form 

A3 4 A3 

/(u, 7 ) =- — ln(sinh(M + 7)) + 



2N 3 ^ e" m7 . , , n . sinh(2m7) 

+ — * sinh 2mu , 

TV z — ' m sinh (07717 ) 

m=l 



— ln(sinh(u + -7)) H -u 

N K 2 3 iV 



2iV 3 * ^ e -m 7 sinh(m 7 ) nn 

+ ^2. Binh(2mix) / ,; (5.20 

J\ z — ' m smh( 07717) 

m=l 



As we can see, the free energy is the sum of the individual contributions of the sites 
in each representation. So, for A3* = (N3 = 0), we obtain again the results of the 
homogeneous case in the representation {3} ({3*}). 

From the free energy, we can obtain the energy density in the fundamental state, 



dS 
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34f 
2 du 



(5.21) 



u=0 



Doing the calculation the result is 



2| N 



coth 7 + - + 4 



+ 



771 = 1 

oo 



coth(| 7 ) + |+4^e- 



-7T17 



777.7 



m=l 



sinh(2m7) 
sinh(3m7) 

sinh(m7) 
sinh(3m7) 



(5.22 ) 



that again is the sum of the individual contributions of each site representations. 

We can apply the results to the alternating case (A3 = N£ = N/2); the free energy 



is 



fialt) ^ 7 ) - in^jnh^ _|_ 7 )) _ 1 l n ( s i n h( w + -7)) 



1 



-m~f 



m=l 



m 



sinh(2mtt 



2 s s 2 
cosh(im7) 



cosh ( I my) 



(5.23 



and the energy density of the fundamental state 



£W = i coth 7 + coth(| 7 ) +|-3^e 



-7717 



777=1 



cosh(^m7) 
cosh §7777) 



(5.24) 



The solutions we have given, have been obtained by taking hyperbolic functions for 
the solutions (|4.2a-d|) and ( |4.7a-(j|) of the YBE. By considering the trigonometric solutions 
of these equations and following the same steps, we find the BA equations, 



sinh(f 



(i) 



*2 ' 



.sinh(t> j: 1 '* + i- 



N 3 



smh(v k + i 



2 - 

(2) , „-7M-^3 



sinh(f^,' 



n 

3=1 

r 

n 





- vf 


- h) 


s 

n 

1=1 


■ Uf (2) 








- 4 ] 


+ h) 


■ Uf (2) 




+ «§) 


■ W (2) 
sinh^ 




- h) 


a 

n 

1=1 


■ ut (2) 
sinh(iY 


-4 2) 


-i}) 


• uf (2) 




+ h) 


• (2) 
sinh(tY 


-4 2 > 


+ <i) 



(5.25a) 
(5.256) 



In this regime, the roots cover all real numbers (—00, 00). Then, defining an analogous 
function 



_ , x . , sinh(a; + io>) 
$(x, a) = i In — 



sinh(a; — ia) ' 

we can solve the problem again by using the Fourier transform 

'+00 ^ sin(A;(§ - a)) 



f dk 

1 r+oo 

PjW = 7T / dke' kX p 3 {k). 
^ J -00 



(5.26) 

(5.27a) 
(5.276) 
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The no-holes hypothesis for the ground state gives us the densities 

N 3 sinh(£; 7 ) N* sinh(A^) 



h{k) ~ ^Vsinh(ffc 7 ) + ~W sinh(ffc 7 )' (5 - 28a) 

iV 3 sinh(/c|) jV 3 * sinh(fc 7 ) 
p2(/eJ " ATsinhd^) + ATsinhd^)' (5 ^ J 



and the free energy becomes 



iV 3 J . f°° rffc sinh(fcn) sinh(fc(f - |)) sinh(fc 7 ) | 

/ (u,7 = —< -lnsm u + 7 +2 / — . ... . J. 3 7 ^ f 

N { Jo k smh(/c|) sinh(fc^) J 

iV 3 * f , , 3 N /- 00 rfA;sinh(/c-u)sinh(/c(f - ^))sinh(/c 7 ) 



+ ^|-.n S in( U+ - 7 ) + 2 y o T ) ^ ) }(5-29) 



The density of energy of the ground state is 



_ 3 (Ns 
2{ N 



- cot 7 + 2 



AT 



cot -7 + 2 



sinh(/c(f - ^))sinh(/c7)' 



sinh(/c 
sinh(/c(-| 



7T 



sinh(fc^ 



2 / ^""V" 2 

^)) sinh(/c7) 



sinh(/c§ ) sinh(/c^) 



(5.30) 



We can specify these magnitudes for the alternating case (N 3 = N£ = N/2); they are 



.(alt), n 1, • / n 1, • / 3 , y 00 dfc smh (feu smh {k § - ^ cosh f 7 
/( = -2 ln -(^)-2lnsm(. + -7) +j / o y sinh(fcf } Lh(^) 

(5.31) 

and 

3, _,3 . 3 r , sinh(/ C (|-i))cosh(^) 



£ a " = _^( cot7 -cot(- 7 )) - - / dfc ^ ^ J-^-. (5.32) 

4 V r K 2 IJJ 2j sinh(fcf)cosh(A;^) 1 J 

We can describe other quantum numbers of the eigenvectors of the transfer matrix. 
Let we define the number operators 

?! =N U - Nu, (5.33a) 

Y 2 =N d - N d , (5.336) 

where 

JV 

N a = ^ !i ® • • • ® 1 ® (w a )i ® 1 <8> • • • ® Ijv, (5-34) 
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and 

"«!">= {I '/ %l- < 5 - 35 > 

The operators Y\ and Yi commute with the transfer matrix 

[%t(6)} = i = l,2. (5.36) 

The commutation relations with the S-operators are 

\Y 1 ,B i (e)]=-B i (e), (5.37a) 

[y 2 ,Bi(<9)] =<^(0) • (5-376) 

Then, if we apply Y\ and Y 2 on the state obtained by the aplication of r operators 

£? to the pseudovacuum state ||0 > in the first step and s operators in the second step, we 
find 

=(N 3 -r)*((t), (5.38a) 

? 2 *(a0 =(r-s)*(/2), (5.386) 

we have the quantum numbers of this problem as 

iV u - jV s =jV 3 - r , (5.39a) 

iV d - N d =r - s, (5.396) 

and obviously 

N s -N- s = s-N;, (5.40) 

being N q the eigenvalues of iV g , (q = u, u, a 1 , d, s, s). 

In the thermodynamic limit the fundamental state is characterized by 

(*)_-/>*-^ 

MAMA = *»+™J , (5.41/,) 



and then 







_N 3 -N* 






3N 


N d - ivy 




_N 3 -N* 






3N 






N 3 - N% 


. N , 




3N 



(5.42a) 
(5.426) 
(5.42c) 
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For A^ 3 = or N£ = we recuperate the no mixing chain results. 
In the alternating chain (N 3 = N£ = N/2) we obtain 

fN u -Nu\ fN d -N rJ \ fN a -N g \ 

rV- 1 )* = (-V*) w = (^-^) N = °- < 543 » 

This proves that the fundamental state is formed by pairs uu or dd or ss. 

These methods can be generalized easily to higher representations of sit (3); the only 
change is in the g functions (that we call source functions): they change according to the 
highest weight of the representation. The generalization to su(n) is straighforward too, 
but laborious; the MBA will have (n — 1) steps and will be described by (n — 1) source 
functions, related by (n — 1) Bethe equations. 
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